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Dynamics of the TrES-2 system 
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The TrES-2 system harbors one planet which was discovered with the transit technique. In this work we investigate the 
dynamical behavior of possible additional, lower-mass planets. We identify the regions where such planets can move on 
stable orbits and show how they depend on the initial eccentricity and inclination. We find, that there are stable regions 
inside and outside the orbit of TrES-2b where additional, smaller planets can move. We also show that those planets can 
have a large orbital inclination which makes a detection with the transit technique very difficult. 
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1 Introduction dynamical properties of the systems with known transiting 

planets. In this work, we investigate the gravitational influ- 
The p lanet TrES-2b was discovered in 2006 dO'Donovan et ali enee of TrES-2b on a potential second planet in the system. 



2006). It orbits a GOV star with approximately the same 
mass and radius as the Sun. The planet has a mass of 1.198 
Mjtjp and a radius of 1.22 Rjijp- Its orbit is circular and 
has a semi-major axis of 0.0367 AU (see table Q]for all pa- 
rameters and the corresponding error bars.) 

TrES-2b was discovered using the transit method where 
one tries to observe the decrease in luminosity of a star 
caused by a transiting planet. This method is very promis- 
ing and current and future space-mission (e.g. CoRoT and 
Kepler) will use this method to largely increase the num- 
ber of known transiting planets. Although it is possible to 
detect quite small planets via the transit technique the dis- 



Depending on the orbital elements of TrES-2b we can iden- 
tify the zones of stability, where it could be possible to de- 
tect an additional planet. 

Our work is motivated by the very fact that for th e TrES-2 

planet , transit time variatio ns h ave been reported b y lO' Donovan et aL^ 
(2006). lRatz etail d2008al) and lRatz et alT d2008b). 



2 Numerical Setup 

The majority of planetary systems consist of a star and a 
giant planet orbiting in an eccentric orbit around it. TrES- 



covery of earth-like planets is still very difficult. The premisesl 2 is such a single-planet system and can be modeled by 



are better in systems were a transiting planet is already dis- 
covered. By investigating the variations in the transit time 
of a planet it is possible to identify also very small plan- 
ets (IFord & Gaudi 20061 iFord&Ho lman 2003). 
It is thus desirable to have a detailed knowledge about the 
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Tabl e 1 Basic data of the TrES-2 system 



from lO'Donovan et alJ (120061) : I Winn et al.l (12008) 
TrES-2b 



Mass 

Semimajor axis 
Orbital period 
Eccentricity 
Radius 

Time of transit 
Inclination 



1.198 ±0.053M JUP 
0.0367 ± 0.0012 AU 
2.470621 ± 0.000017 days 


1.22 ± 0.045i?j UP 
2453957.63479 ± 0.00038 
83.9 ±0.22 



the elliptic restricted three-body problem, a relatively sim- 
ple dynamical model, where two massive bodies, called the 
primaries, move in elliptic orbits around their common cen- 
ter of mass, and a third body of negligible mass moves un- 
der th eir gravitational influence (for details see ISzebehelv 
( 11967b ). For the computation of the stability maps we varied 
both the orbital elements of the giant and the test particles 
and for more than 350 thousands initial conditions the dy- 
namical character of the massless particle representing an 
earth-like planet was determined. To compute the stability 
maps, the method of the maximum eccentricity (ME) and 
the Lyapunov characteristic indicator (LCI) were used as 
main tools for stability investigations of massless terrestrial- 
like planets: 

- To indicate stability, the ME method uses a straightfor- 
ward check based on the eccentricity. This action-like 
variable shows the probability of orbital crossing and 
close encounter of two bodies and therefore its value 
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provides information on the stability of orbits. This sim- 3 Stability Regions 

pie check was already used in several stability investi- 
gations, and was found to be a powerful indicator of the 
stability character of o rbits ( Dvorak et al.l2003 ; Nagy et al 



2006tlSuri et aDl2005h . 



- As a complementary tool, we computed also the LCI, a 
well-known chaos indicator. The LCI is the finite time 
approximation of the largest Lyapu nov exponent (L CE), 
which is described in detail in e.g. iFroeschld (1 1 9 841) and 



Lohin ger et all (1993). 



The mass of /i of the TrES-2b system is approximately 
0.0022097 solar masses which can be rounded to 0.002. 
The time span of the integration was 10 5 revolutions of the 
primaries. For the integration of the system we applied an 
efficient variable-time-step algorithm, known as Bulirsch- 
Stoer integration method. The most important feature of the 
Bulirsch-Stoer algorithm for N-body simulations is that it 
is capable of keeping an upper bound on the local errors 
introduced due to taking finite time-steps by adaptively re- 
ducing the step size when interactions between the particles 
increase in strength. The parameter e which controls the ac- 
curacy of the integration was determined by preliminary test 
runs and was set to 10~ 8 . 

The giant planet's orbital elements are: semi-major axis 
a\, eccentricity e\, argument of periastron u>i, and mean 
anomaly Mj, The semi-major axis of the giant is a\ — 
0.0367 AU, and its eccentricity was 0. Since the giant re- 
volves in a circular orbit, u>\ is undefined and the mean 
anomaly alone determines the position of the giant along 
the orbit. M\ was varied between 0° and 180° with A Mi = 
45°. 

The semi-major axis of the test particles ranges from 
0.014 to 0. 147 AU for stability maps in the a — i plane and 
up to 0.183 AU for stability maps in the a — e plane. The 
inner border for the test particles is very close to the star. Al- 
though earth-mass planets could exist here, due to the large 
temperatures they would certainly not be earth-like planets. 
The step-size Aa = 1.8-10" 4 AU, Ae = 0.02 and Ai = 1°. 
The upper limit of the grid in e is 0.5 and 50° in i. 

To investigate the resonances in detail we also applied 
the restricted three body problem. In this case the integra- 
tion of the Newtonian equations of motion was undertaken 
with a different method. We used the so called Lie int e- 
gration dHanslmeier & Dvorakl 19841 : iLichtene g gerll 1 9 84 )), 
which implements an automatic step-size and is, because of 
the recurrence of the Lie-terms, a very precise integration 
method. Again the integration time was 10 5 revolutions of 
the primaries. For the calculations we took into account the 
resonances up to the 7:4 mean motion resonance for the in- 
ner and the outer region. For the motion in resonances the 
initial starting position is quite important, therefore we have 
chosen 8 starting positions for the fictitious planets in a cer- 
tain resonance with the mean anomaly M ranging from 0° 
to 315° with AM = 45°. 



FigureQ]shows the stability plots for the a — e plane. We in- 
vestigated the small region between the star and the known 
planet and the region outside the orbit of TrES-2b up to 
semimajor-axes of 0.18 AU. The first 5 pictures show the 
stability region for different initial values of the mean ano- 
maly of TrES-2b whereas for the last picture we took the 
global maximum eccentricity that was obtained for any value| 
of Mi. 

The overall shape of the stability is of course not influenced 
by the mean anomaly; it looks the same in all cases. One 
can see that the border of the inner region between the star 
and TrES-2b is almost independent of the initial eccentric- 
ity of the test particle. For initial semimajor axes smaller 
than 0.025 and initial eccentricities smaller than 0.25 the 
test particles stay stable on moderate eccentric orbits with 
e < 0.25. Particles with larger initial eccentricities can still 
stay on stable (eccentric) orbits further away from the planet 
(a < 0.02 AU). 

The border of the stability region outside the orbit of TrES- 
2b shows a different behavior. For almost circular orbits the 
whole range between 0.05 an 0.18 AU shows stable mo- 
tion. This region however is shrinking when the initial ec- 
centricity is increased. The gravitational influence of TrES- 
2b reaches out to 0.12 AU; from this point on the stability 
region is independent from the initial eccentricity. 

The influence of mean-motion resonances with TrES- 
2b is also visible in the figures. Depending on the initial 
value of Mi, one can observe additional stable regions cor- 
responding to the position of the resonances. They are most 
stable for Mi = 90°; an initial value of 45° shows fewer 
stable resonances; for 0°, 135° and 180° almost all reso- 
nances become unstable (for a detailed discussion of the 
stability of resonances see section@]l. 

The orbital inclination of a planet plays an important 
role when one wants to observe a transit. We know that - 
as viewed from the Earth - TrES-2b passes the visible disk 
of the star. Normally one would expect any additional plan- 
ets to orbit almost in the same plane as TrES-2b. But from 
our solar system we know that small objects like asteroids 
can have substantial orbital inclinations and even the dwarf- 
planet Pluto has an orbit that is inclined about 17° towards 
the ecliptic. Knowledge about the possible inclinations that 
additional planets in the TrES-2 system can have is thus cru- 
cial for the search of transits. 

We thus also investigated the stability regions in the a — i 
plane; the results can be seen in figure [2] For different ini- 
tial values of the eccentricity we plotted the maximum ec- 
centricity in the region between 0.014 and 0.147 AU in the 
semimajor axis and 0° and 50° in the inclination. 
Figure [2] shows clearly that the borders of the stability re- 
gions behave almost indepently from the initial inclination. 
For an initial eccentricity of 0.01 the whole region exterior 
to TrES-2b between 0.05 and 0. 147 AU shows ordered mo- 
tion on circular orbits for all initial values of the inclination 
(with the exception of the 1:2 and 1:3 resonance). The in- 
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Fig. 1 Stability plots in the a — e plane for different initial values of M showing the maximum eccentricity (ME). 
From top left to bottom left: M=0°, M=45°, M=90°, M=135°, M=180°. Bottom right: combination of all figures using the 
maximum ME over all values of M. 



terior region is also stable up to 0.03 AU although here the 
orbits become more eccentric for initial inclinations of 30° 
and larger (with the exception of the 2:1 resonance) where 
the test particles can maintain almost circular orbits even 
for larger inclinations. If the initial eccentricity is increased, 
the maximum eccentricity increases for all initial values of 
a and i. The border of the exterior stability region moves 
inward with increasing initial e from 0.05 AU to rj 0.115 
AU for e = 0.5 (again, at some mean motion resonances 
the orbits can remain stable). 

Figure [3] shows, how the borders of the stability region 
depend on the initial eccentricity. The change of the in- 
ner border with the initial eccentricity is almost constant 
whereas the outer border changes linearly (with a slope of 
approximately 0.1). The difference between i — 0° and 



i = 40° is very small and only significant for small initial 
eccentricities. 

4 Stability of Resonances 

Resonances play an important role concerning the possi- 
bility of additional planets. On the one hand, a resonant 
configuration can be more unstable and thus make it pos 



sible to exclude larger planets (see e.g. iMiller-Ricci et al. 



(200833)). On the other hand, resonances can stabilize the 
motion of planets under certain conditions and we know 
about planetar y systems that are in a mean-motion reso- 
nance (see e.g. lFerraz-Mello et al.1 (|2P05) and references therein).| 

We thus decided to perform a detailed study of the sta- 
bility inside the mean motion resonances. The stability plots 
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Fig. 2 Stability plots in the a — i plane for different initial values of e showing the maximum eccentricity (ME). From 
top left to bottom right: e=0.01, e=0.1, e=0.2, e=0.3, e=0.4, e=0.5. 



presented in section|3]give an overview on the stability prop- 
erties - but to decide whether a resonance is really stable or 
not, one has to use special initial conditions. In this simu- 
lation the test particles were placed directly inside the res- 
onance. Since the initial mean anomaly can have an impor- 
tant influence on the stability of a resonance (see figure [TJ 
we performed calculations for 8 initial values of M of the 
test particle between 0° and 315°. 

Figure|4]shows the results of this study. We investigated 
all major inner and outer resonances with the known planet. 
For the outer resonances (figure |4] top) we can see, that all 
resonances between 0.107 AU (1:5) and 0.053 AU (4:7) are 
stable for all values of M. From 0.044 AU (3:4) to 0.041 
(5:6) the resonances are unstable for all values of M. In 
between lie the 3:5 and 2:3 resonance where the stability 
depends on the mean anomaly. 



The inner resonances (figure|4]bottom) show a more ordered 
behavior: from 0.011 AU (6:1) to 0.026 AU (5:3) all reso- 
nances are stable; the resonances close to the planet, from 
0.028 AU (3:2) to 0.032 AU (5:4) are unstable (with the ex- 
ception of the 3:2 resonance for M = 0°). 

5 Conclusions 

We investigated the dynamical stability of the TrES-2 plane- 
tary system in order to identify regions where an additional 
smaller planet could exist. Our results can be summarized 
as follows: 

- There is broad range of possible orbits for additional 
planets on almost circular orbits. Such a planet could 
also exist inside the orbit of TrES-2b. Only the region 
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Fig. 3 Borders of the stability region and their dependence on the initial eccentricity. The borders are shown for i=0° and 
i=40°. 



between approximately 0.03 and 0.05 AU is excluded 
due to the gravitational perturbations of TrES-2b. 

- Additional planets can also exist on moderately eccen- 
tric orbits. Here, the excluded region increases with in- 
creasing initial eccentricity. Planets with eccentric orbits 
(e > 0.4) can only exist at a distance from the star larger 
than 0.095 AU. 

- Concerning the stability, a possible inclination of the or- 
bit of an additional planets plays a less important role 
than its eccentricity. Planets on almost circular orbits 
and a semimajor axis larger than 0.05 AU can have an 
orbital inclination up to 50° without becoming dynam- 
ically unstable; planets inside the orbit of TrES-2b can 
be stable up to 35°. 

- Many of the major mean-motion resonances in the sys- 
tem are stable. Only close to the existing planet, they 
can destabilize the motion of a planet. It is thus possible 
for planets to exist in a resonant configuration. 

According to our results, the possibility to detect an ad- 
ditional planet in the TrES-2 system are good - at least from 
a dynamical point of view. There exists a wide range of re- 
gions where such a planet could exist on a stable orbit. Our 
results also showed that additional planets can have a sub- 
stantial orbital inclination. If this is the case, a search with 
the transit technique would fail, but a detection by the radial 
velocity technique would be possible. In this case precise 
measurements of the transit time of TrES-2b could be used 
to help identify such a planet. 

Future works will be dedicated to determine the limits of the 
masses the additional planets can have. It will then be pos- 
sible to make much better predictions for the observability 
of possible transits. 
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Fig. 4 Stability of Resonances. The figure shows the stability of the resonances for different initial values of M. The 
upper figure shows the exterior resonances with TrES-2b; the lower figure the interior resonances. White points indicate 
stable motion; black points show unstable motion. 



© 2006 WILEY- VCH Verkg GmbH & Co. KGaA, Weinheim 



www.an-journal.org 



